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1. Introduction 

In this work we present a tool to model the ionization caused by cosmic rays (CRs) as a com¬ 
plementary tracer to gamma ray observations. CRs are deflected by interstellar magnetic fields due 
to their charge and thus loose all the information about their origin. To observe cosmic rays in the 
vicinity of their sources one therefore relies on tracers of cosmic rays. The most commonly known 
are gamma rays produced in cosmic ray interactions with matter or radiation fields. Since the most 
dominant species of CRs are protons the sources of cosmic ray protons are the dominant sources of 
cosmic rays. Protons prodcue gamma-rays in inelastic proton proton interactions via the produc¬ 
tion and subsequent decay of neutral pions. Therefore, molecular clouds in the vicinity of cosmic 
ray sources act as barometers for the local cosmic ray flux (dependning on the mass and distance 
to the molecular cloud). In fact observations in high and very high energy gamma rays revealed 
proton acceleration in supernova remnants interacting with molecular clouds |Q, j ! 5|]. In this work 
we investigate a CR tracer in molecular clouds complementary to gamma rays. Here we investigate 
the ionisation rate of CR potons and electrons. In contrast to the gamma ray emission, the ionisa¬ 
tion rate is a single quantity integrated over the whole CR spectrum. Therefore, the ionisation rate 
does not contain direct information on the underlying CR proton spectrum. However, in contrast 
to gamma-ray emission, the ionization rate level is not limited to proton energies above 280 MeV 
being produced mainly by lower energies CRs. Therefore cosmic ray ionisation measurements al¬ 
low to access the low energy end of the cosmic ray proton spectrum in interstellar space, offering a 
unique possibility to measure energies below the pion production threshold. 

Here we present the basic principles that determine the ionisation rate caused by cosmic rays. 
With respect to previous works on this subject, we improved the calculations in two aspects: i) 
we use the relativistic differential cross section and ii) we include a self-consistent calculation of 
ionization due to secondary electrons. The use of differetial cross sections allows a selfconsistent 
calculation of the secondary ionizations. The calculations presented here have been implemented 
in a code which is available online (http://crime.in2p3.fr/) and free available for the scientific com¬ 
munity. The online tool calculates the ionisation rate for a given cosmic ray proton and/or electron 
spectrum and provides some basic informations such as the most ionizing energies of the consid¬ 
ered CR spectrum. 

2. Cross sections, secondary ionisations and propagation 

In the following section we briefly summarize the ionisation cross sections used in our code 
CRIME (Cosmic Ray Interactions In Molecular Environments). In addition we derive a general 
expression describing the secondary ionisations. Figure [I] shows the cross section and secondary 
ionisations for molecular hydrogen. Finally we shortly adress the issue of cosmic ray propagation 
into a molecular cloud. These section describes the building blocks of the CRIME online tool. 

2.1 Electron cross section 

We use the relativistic binary-encounter-dipole model for the differential ionization cross sec¬ 
tion of electrons. The relativistic differental cross section is than given (eq. 19 [[!]] ) as: 


2 



CRIME 


Julian Krause 


/ da 
\dW 


47lala 4 N 


'(ff+Pi+PZW 


1 


Ni/N-2 


t + 1 
+ {2-(Ni/N)} 

1 df 


■ + ■ 


1 


1 + 2 1 ' 


w + l t — wj (l + F/2) 2 

1 1 b' 2 


+ ■ 


■ + ■ 


+ 


L(w+1) 2 ( t-w ) 2 (1 + F/2) 2 J 

Pf 


N{w+l)dw\_ \ 1 — j8 f ‘ 


//2 


-ff-ln(2b') 


(2.1) 


Here w = W /I and t = T/I are the energies in units of the ionization potential where T is the 
energy of the primary electron and W the energy of the secondary electron. t p = T/m, b p = B/m 
and Up = U /m are the primary energy, binding energy (ionization potential) and average orbital 
kinetic energy of an electron in natural units. 
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The dipole oscillator strength ^ depends on the target material and can be approximated as a 
power law series for the materials under consideration. 


df 

dw 


■ + ■ 


■ + ■ 


■ + ■ 


/ 


(2.5) 


(1+w) 3 (1 + w) 4 (1+w) 5 (1 + w) 6 

For Hi the fitting parameters are: c=1.1262, d=6.382, e=-7.8055, and f=2.144. 

Integrating equation |2j] over the secondary electron energy (w) yields the total ionization cross 
section (see also eq. 20 in [|T]]): While the integral is partially analytical, the part of the integral 
involving the dipole oszilator strength can be isolated and is given as: 
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For the cases considered here, where the dipole oszilator strength can be described as a simple 
powerlaw series, the above integral can also be solved analytically. 


2.2 Proton cross section 

For the proton ionization cross section we adopt an analytical expression given by [J3|] . This 
expression is based on the molecular promotion model in the low energy regime and on the classical 
binary encounter approximation for higher energies. 

The model equation for the differential cross section as a function of the ejected electron 
energie (w = E e //) in units of the ionization potential (7) and the reduced proton velocity (v = 
fif) is given as: 
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( 2 . 8 ) 

(2.9) 


The parameter w c determines the point where the exponential suipression of the cross section 
sets in and is in fact composed of two parts, w c = A — B with A = 4u 2 — 2v and B = R/4I. As long 
as we are in the low energy part of the cross section which is described by the molecular 

promotion model. Within the low energy regime w c is independent of the energy of the impacting 
proton and determined by the target material. 

For impact energies A 2> B the cross section is based on the binary encounter model. The 
maximum of the total cross section lies in the range where A r\j B (E ~ 10 5 ). The here used 
definitions of A and B are the ones used in []3[]. It is worth to mention that a more general definition 
of 


A = 4E — jVlE 


( 2 . 10 ) 


with y = 4 is given in the binary encounter theory. However, [|3]] report that y = 2 has been 
found experimentally and is thus the value used in the model equation |2.7| . In the classical regime 
(E = I /2/m’ 2 ) and with energy measured in units of the ionisation potential I, one obtains the 
expression used in equation T9. F\ and F% are described by the following fitting functions: 

Fi=L x +Hi and F 2 = L 2 H 2 /(L 2 + H 2 ) (2.11) 


Hi =Ailn(l +v 2 )/(v 2 + fii/v 2 ), and H 2 = A 2 /v 2 + B 2 /v 4 , (2.12) 

LI =Civ Z)| /(l+£iv Z) ' +4 ), and L2 = C 2 v D2 (2.13) 


The ten parameters A\, ...,E\, A 2 , ...TF and a' are obtaind by fitting the model equation to the 
data. In the case of H 2 as a target they are: A\ = 0.80, B\ = 2.9, Cj = 0.86, D\ = 1.48, E\ = 7.0, 
A 2 = 1.06, B 2 = 4.2, C 2 = 1.39, D 2 = 0.48, and a' = 0.87 [||]. To obtain the total cross section 
equation Ol is numerically integrated. 


Relativistic corrections to the Rudd model The model equation of |[3j] is not correct for rela¬ 
tivistic proton energies. To expand the model to relativistic impact energies we change the velocity 
parameter v to hold in the relativistic limit. First we compare the reduced velocity parameter with 
the classical velocity v c l: 


We find the relation 



(2.14) 
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However, in the proper relativistic treatment the velcovity is given as: 


V re l = yj 1. - (m p c 2 / (E + m p c 2 ) 2 X c 

where c is the speed of light. We therefore substitute the velocity parameter v with 


(2.16) 


v = \Jm p /2 x sjm e /(m p /l) x v rei (2.17) 

to introduce a velocity parameter which holds in the relativistic limit and reduces to the old 


parameter for not relativistic energies. Moreover, we can not use Equation |2.9| as it was derived 
using the classical energy-velocity relation. Analog to the above derivation we find that: 


w' c = (4£-2V7t)//-w2. (2.18) 

These changes give a sufficiently correct approximation of the cross section in the relativistic 
limit. Due to the above changes the total ionization cross section stays constant for highly rela¬ 
tivistic energies as the normalization of the differential cross section depends only on the velocity 
which asymptotically approaches c. In reality the cross section keeps increasing logarithmically in 
energy, introducing a small uncertainty for highly relativistic particles. 


2.3 Secondary ionisations 

Each ionization process produces a free electron(here we neglect the case where both electrons 
are stripped from the Hi molecule). Those secondary electrons can produce further ionizations if 
their kinetic energy is above the ionisation threshold. In this paraghraph we want to determine the 
effect of ionisation from secondary electrons as a function of the primary ionization process. Given 
the differential cross section the production rate of secondary electrons is: 
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(2.19) 


Here n is the number density of target molecules, / is the flux of primary particles and da ^ 
is the differential ionization cross section for the primary particles. In the following prime param¬ 
eters always refer to the quantities of the primary particles. Assuming a steady state solution the 
spectrum of secondary electrons is given as: 


dN s _ dN s 
dEdV ~ dEdVdt 


( 2 . 20 ) 


Where z is the residence time of secondary particles within the region of interest. We shall 
assume that all secondaries do not travel any significant distance from the point of origin before 
loosing all their energy. In this case the residence time corresponds to the energy-loss time, and 
all secondaries exist only at the location where the primary ionization took place. Rewriting the 
energy loss function as L(E), the energy loss per column density, the residence time is given as: 


E E 

1 = Tioss = “7777777 = j ( 77\ (2-21) 

dE/dt L[E)nv 

Where n is the number density of the target material responsible for the energy losses and v the 
velocity of the secondary electrons. For simplicity we consider a molecular cloud consisting of 
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molecular hydrogen only. The energy loss function is taken from |[2j] based on measurements. 
Therefore it accounts for all relevant processes, namely: ionisation losses, losses due the excitation 
electronic levels and vibrational states, as well as Bremsstrahlung. Using the above expression 
for the energy losses and the spectrum of secondaries we obtain the flux of secondary electrons 
assuming isotropy as: 
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( 2 . 22 ) 


The flux of secondary electrons is proportional to the ratio of production rate to loss rate as 
expected in a steady state solution. Now we can calculate the ionization rate for the secondaries as 
we did for the primaries. To better understand the effect of secondaries, it is convenient to define a 
function, <f ) (£' , )> as the average number of secondary ionizations per ionization event of a primary 
with energy E'. 
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We find that < f > (£’ / ) does not depend on the flux of primaries but only on the total and dif¬ 
ferential ionization cross section of the primary species, and the total ionization cross section and 
the energy losses of the secondaries (electrons). Therefore, for a given target material and a given 
energy of primary particle, ( P(E) is a constant. We note that for multi component targets the defi¬ 
nition of <f>(£■') has to be understood as a function of multiple components where <P(E') tp)ts is the 
number of secondary ionizations of target species t s caused by primary ionizations of target species 
t p . Interstellar molecular clouds are heavily dominated by molecular hydrogen which shows a 
higher ionisation cross section than other abundant targets (like Helium). Therefore, secondaries 
produced on any other target than molecular hydrogen will be less in number and mainly lead to 
further ionization of molecular hydrogen. For the same reasons also the energy loss of secondary 
electrons is dominated by molecular hydrogen. As can be seen from the right panel of Figure 1, for 
H2 <f>(£' / ) is of the order of unity for both electrons and protons and thus represents an important 
correction. 


2.4 Propagation of cosmic rays into a molecular cloud 

To estimate the ionisation deep inside a molecular cloud it is important to correctly describe the 
propagation of CRs inside the cloud and to account for their energy losses. While this in general 
a complex field of research we follow [^|] and assume straight line propagation and continuous 
energy losses along the line of sight. This reflects the minimum amount of energy losses to any CR 
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Figure 1: Shown are the total ionisation cross section (left hand side) and the secondary ionisations per 
primary ionisation (right hand side) as a function of the energy for both electrons and protons on a pure Hi 
target. 


spectrum. It is important to note that even though CRs below a certain energy E s can not penetrate 
beyond a certain column density, a low energy tail of the CR spectrum below E s is present at any 
depth. This low energy tail is produced by CRs of energies E s + e and its shape is given by the 
inverse energy loss function (see [||]) for details). The importance of the low energy tail for the 
total ionisation rate depends on the shape of the assumed CR spectrum outside the cloud and the 
penetration depth considered. 


3. CRIME the online tool 

CRIME is a webrun tool accessible under http : / /crime . in2p3 . fr/webrun . html 
that allows the calculation of the CR ionistion rate of Hi. The user can specify an electron and/or 
proton spectrum among a power law or broken power law shape as well as pre-defined mi nu mi m 
and maximum CR background models. In the near future the choice of arbitrary CR spectra will 
be implemented via a file upload (data points). It is also possible to account for the straight line 
propagation of the spectrum to the column density of interest, alternatively a minimum CR energy 
considered can be specified. Based on the user choices, the calculation and the displaying of the 
results takes 10-40 seconds. The output of the calculation is the ionisation rate per second as well as 
two plots showing the user defined CR spectrum (and the propagated spectrum if requested) and the 
differential contribution to the ionisation rate. The region contributing to 90% of the total ionisation 
rate is highlighted. In case the straight line propagation is calculated, also the shielded energy E s 
is calculated and marked in the plot showing the differential contribution to the ionisation rate. 
An example output for a Voyager -like CR proton model and a maximum CR background electron 
model is shown in Figure In this example the ionisation rate per second is 1.8 x 10 17 and 
3.6 x 10- 15 for protons and electrons respectively. This tool will provide a fast and easy access to 
the calculation of the CR ionisation rate as well as basic information like the most relevant energies 
causing it. 
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Figure 2: Example Output of a CRIME run. Left hand side: Shown are the CR input spectra for electrons 
and protons before and after straight line propagation to a column density of Hi of 10 22 cm 2 . The part of the 
spectra marked in red contribute to 90% of the ionisation rate. The right hand side shows the differential 
contribution to the ionisation rate in arbitrary units. The regions contributing to 90% of the ionisation are 
highlighted and in addition the energy E s not able to penetrate to the chosen hh column density of 10 22 cm 2 
is shown. We note that for electrons the low energy tail of the spectrum produced in the cooling process 
accounts for about 50% of the ionisation, while for protons its contribution is minimal. 

4. Conclusion and Outlook 

We provide a easy-to-use web tool to calculate the ionisation rate of CR electrons and pro¬ 
tons based on relativistic cross sections and a selfconsistent calculation of the effect of secondary 
ionisations. Straight line propagation of the CR spectra is available as a mini mum energy loss cal¬ 
culation. In addition to the value of the ionisation rate the tool provides information on the effect 
of the propagation model as well as the differential contribution to the ionisation rate. 

We are currently working on the implementation of the calculation of the gamma-ray emission 
of the specified cosmic ray spectrum (with the pion production channel already available on the 
webpage) to allow a simultaneous calculation of ionisation rate and gamma-ray emission. The 
whole calculation is based on a very modular code which allows an easy implementation of new 
features like more complicated propagation models or other target materials for the ionisation (with 
Helium being already implemented in the underlying code). 
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